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We investigate the time an electronic excitation travels in a supermolecular setup using a mea- 
surement process in an open quantum-system framework. The approach is based on the stochas- 
tic Schrodinger equation and uses a Hamiltonian from time-dependent density functional theory 
(TDDFT). It treats electronic-structure properties and intermolecular coupling on the level of 
TDDFT, while it opens a route to the description of dissipation and relaxation via a bath operator 
that couples to the dipole moment of the density. Within our study, we find that in supermolecular 
setups small deviations of the electronic structure from the perfectly resonant case have only minor 
influence on the pathways of excitation-energy transfer, thus lead to similar transfer times. Yet, 
sizable defects cause notable slowdown of the energy spread. 
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I. INTRODUCTION 

Electronic excitation-energy transfer (EET) after light 
absorption is one of the key processes in the natural light 
harvesting (LH) event and a prerequisite for charge gen- 
eration in the LH reaction center [ll |2] • Hence, it enables 
efficient energy recovery in biological systems. The effi- 
ciency of the light harvesting process is determined by the 
rates of charge transfer and EET of many single trans- 
fer steps contributing to the entire mechanism. In multi- 
chromophoric supermolecules as for instance LH systems, 
such rates are affected by a number of different proper- 
ties and phenomena: the electronic structure of the single 
chromophores, the electronic coupling between different 
system components, the geometry and arrangement of all 
constituents, the energetic and position (dis)order, and 
the interplay with the environment. 

Recent investigations of the natural LH event have 
emphasized the importance of quantum-mechanical ef- 
fects to the LH mechanism PHH] and are exploring the 
existence of long-lived quantum coherences during the 
EET process even at ambient temperature conditions 
[T5] . The character of these coherences, i.e., the question 
whether they are of electronic or vibronic nature, is a 
controversy in the field |17fET| . The quantum-coherent 
energy transfer is understood to direct and hence sup- 
port energy flow within the ring system of LH antenna 
complexes |21 El |9] [12]. Also in artificial 7r-conjugated 
polymers, similar quantum-coherence effects have been 
found |22fE5| . Thus, quantum-mechanical interference 
of energy-transfer pathways plays a prominent role to 
obtain results beyond the classical laws of probability 
[5] and may improve the efficiency of intrachain energy 
migration [22 . Hope is high that a deep understand- 
ing of the fundamental design principles of natural LH 
complexes provides tailored concepts for the design of 
artificial organic devices [T51 E51 [27] . 



The role of coherent energy transfer and the environ- 
ment of the system has been studied theoretically based 
on the density-matrix formalism where usually the mod- 
els of Haken and Stobl, the model of Redfield, polaron 
modifications, or related theories are employed [5] [251 - 
I35| . Such approaches rely on input from electronic- 
structure theory or from the experiment. In the present 
manuscript, we provide an alternative to the previous 
approaches: We address the energy-transfer problem di- 
rectly on the level of electronic-structure theory using 
time-dependent density functional theory (TDDFT) |36f - 
[IT] in real time. The strengths of TDDFT in the context 
of EET are versatile. First, in principle it allows for 
an exact treatment of the electronic structure including 
electronic excitations and couplings between fragments 
of supermolecular setups. Second, TDDFT is able to 
deal with multichromophoric systems of practical rele- 
vance on a first-principles level at bearable computational 
cost. Yet, the validity of TDDFT results strongly de- 
pends on the quality of the underlying approximations 
for exchange-correlation (xc) effects. Therefore, work is 
going on to improve the reliability of currently available 
xc functionals f 41f[15] . Based on a given xc functional, 
TDDFT allows for a transparent analysis of the influ- 
ence of electronic-structure properties on EET in sizable 
molecular systems. 

The standard TDDFT approach rests upon closed 
quantum-system time-dependent Kohn-Sham (TDKS) 
equations that evolve coherently in time. It does not 
comply with experimental situations that involve energy 
dissipation and/or decoherence. Therefore, to investigate 
the real-time dynamics of a system in contact with a dis- 
sipative bath, one needs to go beyond standard TDDFT 
and use open quantum-system schemes [49^645 . The ap- 
proach that we employ in the present work is embedded 
in the framework of TDDFT and uses stochastic TDKS 
equations |50f[55] . 
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The aim of our work is to investigate the coherent EET 
process in a ring configuration of model molecules in the 
spirit of the circular LH complexes. In this setup, we ex- 
ploit the strengths of TDDFT and examine the influence 
of intermolecular coupling as well as energetic disorder 
and mismatch. To this end, we introduce an effective 
bath operator that plays the role of a measurement pro- 
cess to determine the time, a dipolar excitation is trav- 
elling in the molecular assembly before it dissipates. It 
couples to the dipole moment of the supermolecule or 
appropriate subsystems of the latter. We use this mea- 
surement process to deexcite the system, thus remove the 
entire excitation energy. Note, that in real systems exci- 
tation energy transfer may also occur to the environment, 
e.g., to neighboring complexes and/or may dissipate to 
nuclear degrees of freedom. 

The paper is organized as follows. In Sec.|TT]we outline 
the stochastic TDKS approach and introduce a bath op- 
erator that couples to the dipole moment. More insight 
into our single-particle scheme is provided in Appendix 
\K\ The model system investigated in this manuscript is 
introduced in Sec. |III| where we also examine the perfor- 
mance of the bath operator. Moreover, we explain the 
idea of using the bath operator in the sense of a mea- 
surement process to determine the effective time an ex- 
citation can travel in a supermolecular system before de- 
excitation occurs. Details about the electronic structure 
of the subsystems and the intermolecular coupling in the 
supermolecule are compiled in Appendix ^ We intro- 
duce our method to obtain these properties in Appendix 
[C| In Sec. |IV[ we report our results on the influence of 
varying electronic-structure properties and different in- 
termolecular coupling mechanisms on EET times in the 
model system. Finally, we summarize and conclude our 
findings in Sec. \V\ 



II. DISSIPATION IN REAL-TIME 
TIME-DEPENDENT DENSITY FUNCTIONAL 
THEORY 

A. Stochastic Schrodinger equation 

One of the most basic ideas of open quantum-system 
treatments is to split the entire Hamiltonian into system, 
bath, and system-bath interaction contributions 



electronic many-particle Hamiltonian 



H = Hfi <Si In + h (E) Ht^ 
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is and Ib denote identities in the system (S) and bath 
(B) Hilbert spaces. The system includes all dynamics 
and observables of the core system, e.g., of one molecular 
complex. The bath and system-bath coupling describe 
the environment as for example surrounding molecules 
and its interactions with the core system part. Here, 
the system degrees of freedom are coupled to a bosonic 
environment. The system of interest is described by the 



[p, + Aext(rj,i)]' 
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where Acxt is an external vector potential, Uext a scalar 
external potential, and W describes the particle-particle 
interaction. The environment, given by Hbi induces a 
fluctuating force that drives the system due to the inter- 
action 



SB 



E 
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between the system and the bath. The latter is expressed 
by many-particle operators Sa and Ba in the most gen- 
eral case, where Sa operates on the system degrees, 
on the bath degrees of freedom, and a denotes differ- 
ent possible system-bath coupling mechanisms. The bath 
may have a complex structure and not all of its micro- 
scopic details are relevant for the system dynamics. The 
driving force that is induced by the bath may, therefore, 
be subsumed by stochastic fluctuation and dissipation 
contributions that can usually be characterized by mean 
values and correlation functions The parameter A 

determines the strength of the system-bath interaction 
and can be used as an expansion parameter in the weak 
coupling regime. 

The stochastic Schrodinger equation (SSE) and quan- 
tum master equations are approaches to describe the dy- 
namics of such a core quantum system in contact with 
an external bath [5J1 [SSJ . While quantum master 
equation approaches are based on density matrices, the 
SSE uses a statistical ensemble of state vectors to unravel 
the open quantum-system dynamics directly on the level 
of wave functions. A noteworthy difference between be- 
tween the SSE and Redfield and similar master-equation- 
type approaches is related to the statistical operator. In 
the latter methods, positivity of the statistical operator 
may not be guaranteed in cases of stochastic Hamiltoni- 
ans, time-dependent (TD) Hamiltonians, or TD bath op- 
erators fST, 59, 70 . However, the positivity issue does not 
arise in the SSE. It yields normalized wave functions with 
a positive weight for arbitrary TD operators. This aspect 
is particularly important if one intends to use Hamiltoni- 
ans from density functional theory (DFT) |51l [55] as we 
do in the following sections. Therefore, the SSE is a solid 
starting point for a (TD)DFT theory of open quantum 
systems. 

The open quantum-system approach that is applied 
throughout this work is based on the SSE in the so-called 
Born-Markov limit |M1 IS 



(4) 
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where la{t) are stochastic processes with vanishing en- 
semble average and (S-time-correlation 

u^y = 0, ut)W) = 0, W)W) = So.p6it - 1'). (5) 

The bar denotes the statistical average over an ensem- 
ble of stochastic processes. For the sake of convenience, 
the coupling strength parameter A has been absorbed in 
the bath operator Sa ■ The first term of the SSE Q de- 
termines the usual unitary system evolution under the 
action of the system Hamiltonian Hs ■ Although the SSE 
employs the system Hilbert space only, the coupling to 
the bath is still included by the second term that de- 
scribes dissipation effects due to the system-bath inter- 
action. Finally, the third term introduces fluctuations in 
the time evolution: Although the dissipative term causes 
the probability density to decay in time, the norm of the 
state vector ^(t) averaged over a statistical ensemble of 
realizations is conserved up to fourth order in the system- 
bath coupling parameter A. In the following, we focus on 
a single bath operator S and, for the sake of clarity of 
our notation, omit indices at the bath operator. 

As a result of the stochastic nature of Eq. Q, the 
system wave function can not be simulated by a single 
evolution of the SSE but needs to be represented by a sta- 
tistical ensemble of wave functions {^'^(t)}. Accordingly, 
starting from a pure initial state, the full time evolution 
of expectation values 



Os{t) - (Os> = (*(t)|Os|*W) 



(6) 



of physical observables needs to be calculated from the 
statistical average over all ensemble members as indi- 
cated by the bar. A reliable computation of smooth ob- 
servables requires a large enough set of stochastic real- 
izations. 



B. Stochastic Schrodinger equation and 
Kohn-Sham density functional theory 

So far, the complexity of solving the entire system 
Hamiltonian has been reduced via the SSE framework. 
Yet, still the question about the choice of bath opera- 
tors and potential difficulties with describing the core 
system part remain. We address the question of appro- 
priate electronic-structure theory methods for Hg in this 
section before we turn to bath operators in the subse- 
quent section. 

The open quantum-system SSE approach has recently 
been extended [5U1 [5T| to the framework of TD current 
density functional theory (TDCDFT) and its applicabil- 
ity has been discussed |53fl55] . In contrast to standard 
TD(C)DFT, the open quantum-system approach oper- 
ates on ensemble-averaged quantities. It requires the 
ensemble-averaged particle density 



where the density operator is defined as 
n(r) = ^(5(r-f,), 

i 

and the ensemble-averaged current density 



j(r,t) = (j(r,t)), 
where the current operator reads 



(8) 



(9) 



(10) 



and {. , .} denotes the anticommutator bracket. The the- 
orem of stochastic TDCDFT of Refs. [50, 51_ states that 
under reasonable physical conditions for a given and fixed 
bath operator S, many-particle interaction W, and initial 
state 5'(r,t — 0), a one-to-one correspondence between 
the external vector potential Aoxt(r, t) and the ensemble- 
averaged current density j(r, t) exists. Thus, the theorem 
paves the way for a non-interacting KS scheme of open 
quantum-system TDCDFT as it guarantees the unique- 
ness of a KS system that yields the same current density 
as the true interacting system. The corresponding KS 
Hamiltonian reads 



KS 



E 



[pi + Aext(ri,t) -h Axc(fi,t)]^ 



(11) 



with the xc vector potential Axc(r,i) and the Hartree 
potential vii{r,t). The xc vector potential in the open 
quantum-system approach may in general depend on 
j(r,t), the initial state $Ks(r,i = 0), and the bath op- 
erator S. The system dynamics can be calculated from 
the KS Slater determinant $Ks(r,t) evolving according 
to the open system KS equation [50] 

fi9t$KS W = ^KS*KS W - ^^^^$KSW + l{t)S^Ks{t). 

(12) 

For practical calculations however, one needs to rely on 
existing approximations for the xc vector potential as the 
true Axc(r, t), especially in open quantum systems, is not 
known |50| . 

The ensemble-averaged current density - the basic 
property of the stochastic TDCDFT framework - is re- 
lated to the ensemble-averaged density via the continuity 
equation of the open quantum system |51l ITT] 



(13) 



The density modulation that is induced by the bath is 
described by 



n{r,t) = (n(r)), 



(7) 



(14) 
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The proof of the above stated one-to-one correspo ndence 
requires implicitly that J-B(r,f) is a functional of n(r,t) 
and j{r,t) alone and can be solved uniquely to determine 
n(r,t) [531 159j. However, although there is a unique re- 
lation between the density and the current density, the 
proof guarantees only the one-to-one correspondence be- 
tween j(r,t) and the vector potential. The ensemble- 
averaged density of the true and the non-interacting 
system are not necessarily equal. Therefore, the proof 
cannot be extended straight forwardly to a stochastic 
TDDFT framework where a one-to-one correspondence 
of the scalar potential Woxt(r,i) and n(r,t) would be 
needed to set up a non-interacting density-based TDDFT 
KS Hamiltonian. Yet, for applications that are based on 
the time evolution of the density such as the following 
ones, it is of great interest to find an open quantum- 
system approach in the framework of TDDFT. 

Extensions of TDDFT to open quantum systems were 
already performed by Burke et al. [49] and Yuen-Zhou et 
al. [50] based on density matrices and master equations. 
Both strategies aim at the representation of the density 
of the interacting open quantum system by an auxiliary 
system with a different particle-particle interaction. To 
prove the validity of such approaches, one needs to find 
the corresponding external potential of the auxiliary KS 
system that produces the same density under given con- 
ditions. The proofs of Refs. [39] and [SIT are based on 
the reasoning of Refs. |36l 1721 175] and rely on assump- 
tions about the action of the bath operator. Yet, details 
about the consequences of these assumptions, in partic- 
ular with regard to practical calculations and choices of 
the bath operator, are open questions in the field, and 
the use of approximate KS Hamiltonians in master equa- 
tions and related implications have not been fully ex- 
plored. A conclusive extension of the proof of stochas- 
tic TDCDFT that abandons the ensemble-averaged cur- 
rent density and establishes an open quantum-system 
TDDFT scheme based on the ensemble-averaged density 
alone is not available yet. 

Therefore, practical applications of stochastic TDDFT 
so far rely on approximate realizations of the stochas- 
tic formalism in KS TDDFT. In the approach of Refs. 
|52l 153] , the system Hamiltonian is approximated by a 
standard TDDFT Hamiltonian 



^KS - 



+ tiext (f , i ) + UH (l-i , t ) + -Cxc (f i , i ) 



— ■hKs{ii,t) 



(15) 



and used in the open quantum-system KS equation ( 12 ) 
together with a heuristic bath operator |57j. Here, we 
apply an open quantum-system KS scheme of the same 
kind and treat xc effects on the level of the adiabatic 
local density approximation. More background about the 
difficulty of establishing a DFT Hamiltonian in the open 
quantum-system framework, as for instance about issues 
like the u-representability problem of the current density 



|74| . is provided elsewhere |75| . 



C. Bath operators in stochastic time-dependent 
density functional theory 

The choice of the bath operator is the second ingredient 
for defining the open quantum-system scheme. Here, we 
apply a bath operator that introduces dissipation of the 
excited quantum system back to its ground state. The 
operator is a projector that reads 



S = ^M[n{r,t)]\<^Ks{to)){<fKs{t)\. 



(16) 



It exhibits the scaling factors ^ and M[n{r,t)]. Choos- 
ing M = 1 amounts to a deexcitation of the entire sys- 
tem with decay rate 7 and decay time constant r = I/7. 
The bath operator ( 16 1 covers a broader range of physi- 



cal applications if one chooses a density-dependent factor 
M[n(r, t)] that renders the operator sensitive to moments 
of the density. For physical processes, as for instance 
absorption and emission, it is interesting to consider a 
dissipation mechanism that depends on the (transition) 
dipole moment of the excited system. Such a coupling 
to the dipole moment can be obtained by choosing M 
proportional to the time variation of the dipole moment 
d{t) = / rn{r,t)d^r, i.e., M[n{r,t)] cx \d{t) - d(io)|. In 
the same way, a coupling to other moments of the den- 
sity, as for instance particle-number, dipole-moment, or 
quadrupole-moment fluctuations, can be introduced via 
M. Note, that we model only the system part of the 
system-bath interaction as the bath side of the coupling 
mechanism is subsumed in the effective decay rate. 



D. Simulation algorithm 

Having introduced the stochastic TDDFT formalism 
and having discussed physically interesting bath opera- 
tors, we finally turn our attention to the practical algo- 
rithm that we used to solve the open system KS equa- 
tions. Our simulations are conducted with the quantum- 
jump algorithm [65, 67, 76, 77, that has been introduced 
in the context of open quantum-system KS equations in 
Ref . [53j . The algorithm is implemented in the Bayreuth 
version |78|[75] of the PARSEC program package ^80^ and 
rests upon the single-particle bath operator 



= ^M[n{r,t)]\^,{to)){^^{t)\. 



(17) 



We demonstrate in Appendix [A]that using this bath oper- 
ator in the following set of quantum-jump single-particle 
equations is equivalent to using S in the corresponding 
algorithm for solving Eq. (12 1. The quantum-jump algo- 



rithm relies on a piecewise deterministic evolution |67| of 
the set of norm-preserving Schrodinger equations 



(18) 
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that are interrupted by quantum jumps. The quantum 
jumps represent the non-deterministic action of the bath 
operator. The points of time where such a jump oc- 
curs are determined by a random process according to 
a waiting-time distribution. However, as this waiting- 
time distribution is not known beforehand, it needs to 
be determined alongside the actual propagation of the 
open quantum KS system (for details, see Ref. The 
waiting-time distribution can be determined from the de- 
cay of the norm 



|^r(r,i)pdV 



(19) 



of an auxiliary system of N particles in contact with the 
bath evolving according to 



'it) - l^u.^r 



it). 



(20) 



In this system, norm conservation is not build in explic- 
itly. One obtains single waiting times by drawing ran- 
dom numbers in the interval [0,1] and choosing the time 
T when r]{T) drops below this number as the quantum- 
jump time. The waiting-time distribution can be calcu- 
lated from many samples of such single jump times. 

Due to the special choice of the bath operator that we 
consider in this work, some practical advantages occur in 
the algorithm: As the bath mechanism relaxes the system 
to its ground state and if no additional external perturba- 
tions act after the initial excitation, the system remains 
in the ground state after each quantum jump. As a con- 
sequence of this model bath interaction, the time evolu- 
tions of all orbitals of the stochastic ensemble follow the 
same pattern: The KS system evolves deterministically 
until a quantum jump occurs. Then, back in the ground 
state, the system propagates trivially with TD phases ac- 
cording to the KS eigenvalues as the bath operator does 
not act on the ground state. Therefore, the full statisti- 
cal ensemble can be generated from a single deterministic 
evolution of Eq. ( 18 ) alongside with the determination of 
the waiting-time distribution: (i) calculate a long enough 
deterministic evolution together with the norm decay of 
the auxiliary system, (ii) draw a large number of random 
numbers, find the corresponding quantum-jump times, 
and determine the waiting-time distribution, (iii) gener- 
ate single trajectories of the ensemble where the time evo- 
lution up to the jump is determined by the deterministic 
evolution and ground-state values are used for the time 
after the jump. The physical quantities of interest follow 
from averaging over corresponding observables according 
to Eq. (|6| calculated from the thus obtained statistical 
ensemble. Moreover, the bath operator preserves the or- 
thogonality of the KS orbitals and, therefore, an explicit 
orthogonalization of the orbitals as it is performed in Ref. 
[53j is not needed in the present scheme. 
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FIG. 1. Illustration of the model setup which is a circular 
arrangement of eight molecules ml, m2, . . . , m8 with equal 
distances between all molecules. The excitation of the entire 
ring setup is performed via a localized momentum boost in 
the region of molecule m5. We measure the excitation spread 
in this arrangement by application of a dissipative bath op- 
erator in terms of a measurement process at molecule ml. 
For a transparent investigation of the influence of the in- 
termolecular coupling and the energetic arrangement of the 
molecules on energy-transfer pathways, we introduce defects 
in the molecules m3 and m7. All other molecules are fixed 
according to their ground-state geometry, but their electronic 
state is allowed to evolve in time. (See App end ix [b| for more 
details about the specific setup used in Sec. IV ) 



III. 



MODEL SYSTEM AND BATH COUPLING 
MECHANISM 

A. Model setup 



Having established the conceptual details of our ap- 
proach, we turn to practical applications. In the follow- 
ing, we deal with a model system of circularly arranged 
molecules that is designed in the style of circular LH com- 
plexes of the antenna system of LH organisms (see Fig.jlJ. 
Our aim is to study the influence of electronic-structure 
properties on the EET time scale and pathways in such 
supermolecules. Here, we focus on the influence of the 
coupling and the energetic arrangement. With this in 
mind, the physical picture behind our calculations is the 
following: One molecule of the complex is excited due to 
light absorption. Then, the excitation travels in the sys- 
tem where the transfer mechanism is coherent EET via 
the electronic degrees of freedom. As a result the entire 
complex gets excited and we intend to measure the time 
scale of this excitation spread. 

In our basic model setup consisting of eight molecules 
as indicated in Fig. [T] all molecules in the circular ar- 
rangement and all intermolecular distances are equal. 
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Initially, starting from the ground state of the setup we 
introduce an excitation at one of the molecules. To be 
explicit in our assignment, we take molecule m5. We 
simulate the excitation due to light absorption by a mo- 
mentum boost. In practical calculations, this boost is 
applied only in one section of the real-space grid, so that 
only molecule m5 gets excited. 

The bath operator plays the role of a measurement pro- 
cess that determines the time an initial excitation trav- 
els within the ring system. In our setup, it operates on 
molecule ml on the opposite side of the initial excita- 
tion to measure the time scale the excitation needs to 
travel halfway through the ring. The measurement pro- 
cess is supposed to model deexcitation, thus remove the 
entire excitation energy out of the system when the initial 
excitation has reached molecule ml. In contrast to true 
molecules where typically due to relaxation the excitation 
moves to lower lying energy levels, here the measurement 
processes immediately brings the supermolecule back to 
its ground state. It is an incoherent mechanism without 
back transfer from the bath to the system. We assume 
that the bath is sensitive to dipolar excitations and cou- 
ples to the dipole moment of molecule ml. For these 
reasons, we take the bath mechanism introduced in Eq. 



( 16 1 in its single-particle version of Eq. ( 17 1 and design it 



to be sensitive to dipolar excitations: It detects changes 
of the dipole moment |di(t) — di(io)| of molecule ml only. 
Here, the index denotes molecule ml and indicates that 
the dipole moment is calculated only in the correspond- 
ing section of the grid. The corresponding scaling factor 
M of the bath operator includes a normalization factor 
D and reads 



M[n{v,t)] 



|di(t)-difa)| 
D 



(21) 



In this approach, the time constant related to the rate of 
the dissipative mechanism is a free parameter, whereas D 
needs to be chosen reasonably as discussed in Sec. |IIIB[ 
In summary, the bath mechanism measures locally how 
much of the excitation energy has reached molecule ml 
and deexcites the entire system accordingly. 

For a transparent analysis of the influence of the 
coupling mechanism and the energetic arrangement, we 
choose dimers to be the components of the circular sys- 
tem (see Appendix [B|. In dimers, the electronic structure 
can be modifled easily by bond-length variation. Thus, it 
is straight forward to introduce defects in terms of mod- 
iflcations of the electronic structure as indicated in Fig. 
[T] To guarantee uniqueness and transparency of our in- 
vestigations in Sec. |IV[ we modify only molecules m3 and 
m7, and flx all other system components. The electronic- 
structure properties of the model system are discussed in 
detail in Appendix [B] However, the following considera- 
tions are not restricted to dimers but can be applied to 
more general systems. In a circular arrangement of model 
molecules, there are at least three time scales relevant for 
the signal observed on the ring: the time scale related to 
the energy of the excited states, the time scale due to 
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FIG. 2. Norm decay 'q{t) of a single model molecule. The 
damping is performed with different decay-time constants r 
after an initial boost excitation with 0.001 eV excitation en- 
ergy. The norm decay always follows an exponential function 
exp(— t/r) with the preset time t. Here, quantum jumps were 
performed as the norm dropped below 0.014 %. They mani- 
fest in vertical lines where the norm jumps back to one and 
stays for the rest of the time evolution as the bath operator 
does not act on the ground state. 



intermolecular coupling (here, more time scales are in- 
volved if the system is partly or in total off-resonant), 
and the time scale due to dissipative bath action. The 
influence of the bath mechanism on a single dimer is dis- 
cussed in the following before we present our results on 
the life time of the excitation on the ring in Sec. |IV[ 



B. Assessment of the bath operator 

In this section, we discuss the behaviour of our model 
bath operator, i.e., we adjust the normalization factor D 



of Eq. (21 1 and investigate the effect of the bath operator 



on a single molecule. The aim of this investigation is to 
find D such that the functionality of the scaling factors 7 
and (|di(i)— di(io)|)/-C' in front of the projector is clearly 
split into two contributions: D needs to be determined 
such that the decay time r of the relaxation process is 
determined only by the decay rate 7, thus t = I/7. The 
dipole dependent scaling factor is intended to account 
for the coupling to the dipole moment of molecule ml 
and should not interfere with the role of 7. Therefore, 
D needs to be adapted to the dipole oscillations of the 
isolated model molecule ml. In the following, we outline 
a procedure for the determination of D and investigate 
some relevant influencing factors. 

To determine D, we calculated 100 fs of the dipole- 
moment time evolution of a single molecule in a closed 
quantum system after an initial momentum boost. We 
apply the momentum boost along the bond axis of the 
dimer. Thus, only the dipole moment along this axis 
gets excited and the investigation can be restricted to 
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20 30 
time [fs] 

FIG. 3. Ensemble-averaged dipole moment (^-component) 
and total energy of a single molecule in contact with a dissi- 
pative bath that induces a decay time of 10 fs. 



this dipole- moment component. There are various differ- 
ent options to compute D based on the existing dipole- 
moment data: It may be chosen to be e.g. the first maxi- 
mum, the average over all maxima, the absolute average, 
or the absolute square average of the dipole moment. We 
find by numerical tests that only normalization accord- 
ing to the average of the absolute square of the dipole 
moment gives an exponential decay with time constant r 
that we aimed at by setting the decay rate. Moreover, we 
checked this observation with different decay-time con- 
stants in Fig. |2] where we plotted the norm decay of the 
auxiliary system of the quantum-jump algorithm. All 
computed decay times complied with the preset decay 
rate. 

Another influencing factor on the determination of D 
is the energy that we introduce into the system by the 
boost excitation. We interpret the excitation procedure 
as a single absorption process of an energy portion. The 
dissipation of such an energy portion due to the measure- 
ment process should be independent of the boost strength 
as long as the boost strength corresponds to a single ex- 
citation process. Therefore, in the single excitation case, 
the normalization factor needs to be adapted to the boost 
strength. We guarantee for the boost strength adaptation 
by always determining D from a closed quantum-system 
calculation with the same initial boost excitation as in 
the open quantum-system calculation we aim at. We 
find that the thus obtained normalization factor yields a 
decay time that is independent of the boost strength. 

The excitation process could likewise be performed by 
application of an external laser field. In this case, an 
instantaneous excitation could be simulated by a short 
pulse and the determination of the normalization factor 
should then be performed using the dipole signal after 
the pulse excitation. 

Having investigated the norm decay in detail, we fi- 
nally give some results on other important observables. 
The exponential decay of the norm translates into an 
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time [fs] 
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FIG. 4. Coherent time evolution of the dipole moment (2- 
component) of four molecules of a circular setup of eight iden- 
tical molecules arranged according to Fig. [l] The intermolec- 
ular distance is 20 bohr. Here, an initial boost excitation was 
applied at f = in the subsystem m5 and no dissipative mech- 
anism acts within the entire system. The oscillation pattern 
emerges due to interference of the dipolar excitation that is 
traveling in the ring system along both directions. 



exponential decay of the total energy and the envelope 
of the dipole-moment oscillation (see Fig. [3|. Thus, the 
bath operator fulfills all desired criteria and is ready for 
application in more complex molecular systems. 



IV. RESULTS AND DISCUSSION 

A. Resonant excitation spread and decay time 
constants 

As a first step towards the investigation of EET in 
circular supermolecules, we study the perfectly resonant 
coupling situation where all molecules are the same. We 
start with the fully coherent case without system-bath 
coupling. Apart from this, we take the setup of Fig. 
[T] with an intermolecular distance of 20 bohr and per- 
form the boost excitation at molecule m5 as explained 
in Sec. |III[ The time evolution of the dipole moment of 
four molecules of this ring configuration is depicted in 
Fig. |4j We observe a fast oscillation of the dipole mo- 
ment that corresponds to the lowest excitation energy 
of our model system at 2.1 eV while the envelope fol- 
lows an interference pattern: At different points of time 
the largest dipole moment amplitude can be observed at 
different subsystems. This interference pattern emerges 
as the dipolar excitation travels in the ring system along 
both directions. It is governed by the intermolecular cou- 
pling strength, i.e., the pattern is determined by time an 
excitation needs to be transferred between neighboring 
molecules. In our model system at a separation of 20 
bohr the coupling strength is 0.038 eV fSTl. If two of 
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FIG. 5. Ensemble-averaged dipole moment (z-component) 
time evolution of four molecules that are arranged as in Fig. 
|4] After the initial boost excitation at molecule m5, the exci- 
tation travels in the ring and decays due to a dissipative bath 
that is acting with a decay time of 5 fs on molecule ml. 



those model molecules were isolated, this would amount 
to a resonance oscillation with a cycle duration of 107.8 
fs. In this case, after each quarter of this period the max- 
imum of the dipole oscillation can be observed on one of 
the two neighbors [81j. Yet back in the circular setup, 
each molecule has neighbors on both sides, thus provides 
a pathway for excitation spread. 

We now add the bath mechanism of Sec. |III| acting on 
molecule ml. To measure the EET time, the bath needs 
to break the interference of the so far fully coherent en- 
ergy spread, thus it needs to operate on a comparably 
short time scale. Here, it acts with a decay time of 5 
fs. This choice will be motivated further in the follow- 
ing. Again, we depict in Fig. [5] the time evolution of 
the dipole moment of the same four molecules. Due to 
the bath mechanism, an additional time scale comes into 
play and one observes clearly how the oscillation of the 
dipole moment decays in all subsystems of the ring due to 
the dissipative process. The coupling to the bath breaks 
the coherent evolution and causes relaxation of the sys- 
tem back to its ground state. To be able to measure the 
traveling time of an excitation in the circular setup, the 
decay time needs to be chosen short enough to sample 
the initial stage of the coherent EET before the inter- 
ference pattern starts to build up. In the example of 
Fig. |5] the decay time of 5 fs fulfills this criterion: The 
dipole oscillation at molecule ml reaches only one maxi- 
mum and subsequently decays to the ground state, thus 
no interference emerges. 

In this open quantum-system simulation, the decay 
of the norm r](T) of the auxiliary system indicates how 
much of the excitation energy has reached molecule ml . 
Therefore, it serves as a measure for the life time of 
the excitation on the ring traveling from molecule m5 
to molecule ml before relaxation occurs. We monitor 
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mi 




rj = U.IUU 


rj = u.UoO 


77 = 0. 


0.0 bohr 


0.0 


bohr 


52 fs 


75 fs 


117 fs 


-0.1 bohr 


0.0 


bohr 


56 fs 


110 fs 


121 fs 


-0.5 bohr 


0.0 


bohr 


84 fs 


97 fs 


143 fs 


-0.1 bohr 


-0.1 


bohr 


57 fs 


112 fs 


123 fs 


-0.1 bohr 


-0.5 


bohr 


85 fs 


97 fs 


157 fs 


-0.5 bohr 


-0.5 


bohr 


166 fs 


181 fs 


211 fs 


removed 


0.0 


bohr 


80 fs 


89 fs 


125 fs 


removed 


-0.1 


bohr 


99 fs 


125 fs 


147 fs 


removed 


-0.5 


bohr 


205 fs 


227 fs 


265 fs 



TABLE I. Times T where the norm ri{T) of the decaying aux- 
iliary system drops below 0.100, 0.050, and 0.012 for different 
modifications of the molecular arrangement. The calculations 
are performed in a circular setup of model molecules (see Fig. 
[T]) with an intermolecular distance of 20 bohr. We modified 
the electronic structure by systematically introducing defects 
in the molecules m3 and m7 in terms of bond-length variation 
and/or removing one of the molecules. These modifications 
are indicated in the first two columns of the table. In our 
model system, a bond-length reduction of 0.1 bohr (0.5 bohr) 
amounts to an energetic detuning of the excitation energy by 
0.049 eV (0.135 eV) and a reduction of the coupling-matrix 
element by 0.006 eV (0.012 eV). 



ri{T) and take the times T when ri{T) drops below pre- 
set thresholds to quantify this life time. However, due to 
the interference along the excitation spread, monitoring 
just a single threshold of this norm decay may be domi- 
nated by interference effects, thus its meaningfulness may 
be compromised. Therefore, we considered the norm de- 
cay down to three different thresholds: 0.100, 0.050, and 
0.012. The times T for the perfectly resonant situation 
are listed in Table [T] They serve as a starting point for 
investigating the influence of electronic-structure prop- 
erties on the time an excitation needs to spread in the 
circular arrangement of the model system with an inter- 
molecular distance of 20 bohr. 



B. Influence of energetic disorder 

Having examined the resonant excitation-energy 
spread, we now assess the influence of energetic disor- 
der on the EET time scales by introducing defects in 
the ring setup. Such defects can be inserted via bond- 
length variations and/or removing one of the molecules. 
We compare in Table |T] the times T to the fully reso- 
nant situation of the previous section. We find that small 
changes of the bond length of -0.1 bohr in molecule m3 
or in both molecules m3 and m7 result in a minor in- 
crease of the time scale set by the resonant case. Larger 
increases of the time T can be observed in cases where 
in one of the molecules the bond length was changed by 
-0.5 bohr. Note, that already in the data set discussed so 
far the infiuence of interference effects on the life times 
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of the excitation on the ring can be observed: While the 
times corresponding to the decay down to 77 = 0.100 and 
77 = 0.012 give a clear picture, the times T aX r] — 0.050 
exhibit a different behavior. Yet, we observe a clear rise 
of the transfer time as soon as the more severe defect of 
-0.5 bohr bond- length variation was introduced in both 
pathways between the initially excited molecule and the 
molecule where the measurement process works. The in- 
crease of this life time of the excitation on the ring is on 
the order of a factor of two. 

As another defect, we investigated a setup where 
molecule m3 was removed from the ring. In this case, 
the coupling between molecules m2 and m4 is smaller by 
about an order of magnitude than the previous coupling 
strength between m3 and its neighbors [8T . In this situa- 
tion, as long as all other molecules are perfectly resonant, 
the times T increase due to less severe defects in molecule 
m3 by about the same amount as before. Introducing ad- 
ditional defects in molecule m7, the energy-transfer pro- 
cess slows down notably as we observe the longest times 
T in case of removed molecule m3 and bond-length vari- 
ation of -0.5 bohr in molecule m7. 

In conclusion, we find that in cases of resonance or close 
to resonance the EET time scale is not affected much 
by variations of the coupling strength and interference 
effects play a role for the life time of the excitation on the 
ring corresponding to specific setups. However, as soon as 
sizable defects occur in both pathways of energy transfer, 
the time scales determined by our measurement process 
notably increase, thus EET is hindered noticeably. 



C. Influence of the intrasystem coupling 

In the final investigation, we examine the influence of 
the intermolecular electronic-coupling mechanism on the 
life time of the excitation on the ring. Here, we focus 
on the difference between full electronic coupling and the 
dipole-dipole coupling that is performed in the promi- 
nent Forster theory of EET |81fl83| . In order to observe 
noticeable differences between Forster-type and full cou- 
pling, we had to go to an intermolecular distance regime 
where such differences occur. For our model system, it is 
known that significant differences from the dipole-dipole 
coupling can be observed for distances below 20 bohr 
where the dipole-dipole coupling approximation overes- 
timates the true coupling strength ||81J- Therefore, we 
switched to a circular setup of 12 bohr distance between 
the molecules to guarantee for a clear analysis of the in- 
fluence of the coupling behavior. We decided to keep 
the dissipation time constant of 5 fs of our previous cal- 
culations although for 12 bohr intermolecular distance 
the coupling-matrix element of 0.105 eV is notably larger 
than in the previous case [ST]. Our results are compiled 
in Table n 

We find that Forster-type coupling overestimates the 
times T determined by the measurement process in com- 
parison to the true electronic coupling, thus it overes- 



coupling m3 77 = 0.100 = 0.050 77 = 0.012 



full 


0.0 bohr 


95 fs 


116 fs 


163 fs 


Forster 


0.0 bohr 


86 fs 


104 fs 


152 fs 


full 


-0.1 bohr 


99 fs 


140 fs 


182 fs 


Forster 


-0.1 bohr 


94 fs 


108 fs 


154 fs 


full 


-0.5 bohr 


109 fs 


140 fs 


209 fs 


Forster 


-0.5 bohr 


89 fs 


115 fs 


185 fs 



TABLE 11. Investigation of the times T where 77(r) drops be- 
low the thresholds of 0.100, 0.050, and 0.012. The molecules 
are arranged at a distance of 12 bohr. Thus, notable differ- 
ences of the coupling strength exist between the cases of full 
electronic coupling and Forster-type dipole-dipole coupling 
[81| . Moreover, defects in terms of bond-length variation are 
introduced in molecule m3 as indicated in the second column. 



timates EET in the setup. This finding parallels the 
fact that also the coupling-matrix element is overesti- 
mated by the dipole-dipole coupling [8IJ . This result was 
consistently confirmed also in situations where defects 
were introduced in the circular arrangement. Finally, we 
note that the times T observed for an intermolecular dis- 
tance of 12 bohr are smaller than the numbers for 20 
bohr, although in the latter case the coupling is notably 
weaker. This finding comes as a result of interference 
effects that are more pronounced in cases of stronger 
coupling. Therefore, it is difficult to compare absolute 
time scales unless the dissipation time constant of the 
measurement process is chosen short enough to prevent 
a buildup of interference due to excitation-energy spread 
along different available pathways. 



V. SUMMARY AND CONCLUSIONS 

In summary, we used an open quantum-system scheme 
in the framework of TDDFT to investigate the excitation- 
energy spread in a circular arrangement of molecules in 
the spirit of LH complexes. To this end, we introduced 
a bath mechanism that couples to the dipole moment of 
speciflc subsystems of the supermolecular complex and 
operates in terms of a measurement process. We flnd 
that this dissipative mechanism breaks the interference 
pattern of coherent EET and removes excitation energy 
from the system. The time scale of this relaxation is 
determined by the intermolecular coupling and the dissi- 
pation time constant. In this context, it is the strength of 
(TD)DFT to be accurate on the level of available xc func- 
tional approximations in the determination of the elec- 
tronic structure of the entire complex and its subsystems. 
Therefore, the influence of electronic-structure properties 
and intermolecular coupling mechanisms on the EET pro- 
cess can be investigated straight forwardly. Our results 
demonstrate that small deviations of the electronic struc- 
ture, e.g. via bond-length variation, cause only minor 
changes of the time scale monitored by the measurement 
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process. Only when all energy-transfer pathways in the 
circular setup suffer from sizable defects, the EET time 
drops notably. Similarly, variations of the coupling mech- 
anism, as for example via Forster-type dipole-dipole cou- 
pling, change the EET times according to the changes of 
the coupling strength. 
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Appendix A: Single-particle representation of the 
bath operator 

In this appendix, we demonstrate that using the single- 
particle operator Si of Eq. ( 17 ) in the set of single-particle 
equations of Sec. II D is equivalent to solving Eq. (12) 



with the bath operator S of Eq. ( 16 1 using the quantum- 
jump algorithm. In the latter case, one needs to propa- 
gate simultaneously the norm-preserving equation 

fi9t$KS = i?KS$KS - ^^^5$KS + ^||^*KS||'*KS (Al) 



and the auxiliary equation 



aux 
KS 



i?KS<I" 



KS'i'KS 



aux 
KS ■ 



(A2) 



Inserting S of Eq. (16 1 in Eq. (All yields the closed 



quantum-system KS equation 



(A3) 



The latter equation can be translated to the well-known 
set of single-particle KS equations that one also obtains 
by inserting of Eq. (17) into Eq. ([T8|. 



Therefore, it remains to be checked if, given our choice 



of bath operators, Eq. ( A2 1 translated to a set of single- 
particle equations coincides with Eq. (20). To this end, 
we use the method of Ref . |S3] to find the single-particle 



form of Eq. (A2|. The method starts from the reduced 



one-body density matrix 



Pir,r',t) = J $Ks(r,r2, . . . ,rAr) 

$Ks(r',r2,.--,r7v)dV2...dVAr. 



(A4) 



Computing the equation of motion of the one-body den- 
sity matrix by inserting the open quantum -svs tem KS 
equation and the specific choice of S of Eq. ( 16 1 yields 



fiStp(r,r',i) 
/iKs(r,t) - 



p(r,r',i). 
(A5) 



Note, that M is defined in the context of Eq. ( 16 1. Next 



we express p{r,r',t) in terms of a representation of N 
orthonormal basis functions ifj [SI] according to 



JV 



p(r,r',t)=^^,(r,i)^*(r',t). 



(A6) 



Applying this representation to Eq. (A5) yields for the 
case without dissipation (7 = 0) 



N 



N ^ 



(A7) 



For the limit 7 = 0, a single-particle representation can 
be obtained by introducing the time-dependent matrix 



(A8) 



Equation (A7l demonstrates that the matrix Cjj'{t) is 
Hermitian (see also Ref. [53]) and, therefore, can be diag- 
onalized. Based on this finding, we assume that {cpj (r, t)} 
is the diagonal basis, without loss of generality, and that 
it fulfills the eigenvalue equation 



Qdt - hKs{r,t) ipj{t) = ej{t)(fij{t) 



(A9) 



with real- valued eigenvalues ej{t). Next, we consider the 
case which includes dissipation. For 7 7^ we can intro- 
duce in a similar way a matrix Djji(t) according to 



^Ks(r,i)- ^7M^ 



(AlO) 



which includes the dissipative term. Note that Djji (t) 
contains the Hermitian contribution Cjji (t) and an anti- 
Hermitian part ^^AI'^ {(fj{t)\(fji{t)) . Since our dissipa- 
tion operator is proportional to a unit operator, the ma- 
trix Djji(t) also remains diagonal in the basis {<^j(r, t)}. 
While the eigenvalues ej (t) of Cjj' (t) have been all real- 
valued, the eigenvalues ej{t) of Djji(t) aquire a constant 
shift in the imaginary part 



Finally, we can set 

if>j{t) = exp <! -fi 



6,(r)dr^^— (0 



to obtain the set of single-particle equations 

fi 
2 



(All) 



(A12) 



(A13) 
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FIG. 6. Upper panel: Excitation-energy variation (energetic 
olf-resonance) of one Na2 depending on variations of the bond 
length. We computed this property from a supersystem cal- 
culation as explained in Appendix[C]using two different dimer 
distances. Lower panel: Coupling-matrix element calculated 
in a supersystem of two Na2 depending on the energetic off- 
resonance of one of the dimers. 



Note that due to the dissipation, the orbitals ip^^^it) stay 
not normalized during the time-evolution. 
To conclude, we have illustrated that the single particle 
equations in (A12) are equivalent to the set of equations 

is inserted. Therefore, the 



(20 1 when the operator Sj is inserted. Therefore 
validity of the single-particle scheme for solving Eq. ( 12 ) 
with the bath operator S is established. 



Appendix B: Details of the model system 

To guarantee for a transparent analysis of the coupling 
and energetic arrangement, we choose sodium dimers as a 
well-approved model molecule |81l [55] . It exhibits strong 
dipolar character and the electronic structure can be 
modified easily by bond-length variation starting from 
the experimental bond length of 5.78 bohr. All sodium 
dimers are aligned along the z-axis according to the setup 
of Fig. [T| and their centers of mass are in the x-y-plane. 
In the following, we introduce the electronic-structure 
properties of our model system and investigate the inter- 
molecular coupling situation. 



energy of the one prominent peak with variations of the 
bond length can be calculated along two different ways: 
either from straight forward TDDFT propagation of a 
dimer with shifted bond length or in a supersystem setup 
with two dimers where one of the dimer bond lengths is 
shifted (see Appendix |C| . We performed both methods 
and found reasonable agreement of the results shown in 
the upper panel of Fig. [6j We find that due to shifts of 
the bond length on the order of 1 bohr, excitation-energy 
fluctuations of about 0.25 eV can be induced. 



2. Resonant and off-resonant coupling 

Second, we consider the coupling between two neigh- 
boring molecules and therefore isolate a system of two 
Na2. Here, we address two influencing factors, namely 
the distance between the dimers and the energetic ar- 
rangement determined by the dimer bond length. The 
distance dependence of the coupling-matrix element in 
a resonant situation with sodium dimers of equal bond 
length was already investigated in an earlier work [BlJ. 
This investigation was based on the Davydov splitting 
in a resonant two-level interpretation extracted from the 
dipole moment time evolution: For distances above 25 
bohr, the coupling is of dipole-dipole type, whereas clear 
deviations from the dipole-dipole character can be ob- 
served for smaller distances. 

In an off-resonant coupling situation between isolated 
excited states, a similar two-level picture can be 
applied to determine the coupling-matrix element. In 
analogy to the resonant case, also in this situation, 
the coupling strength can be extracted from a suitable 
observation of the dipole moment of separate system 
parts (for details, see Appendix [C| . The results of the 
coupling-matrix element as a function of the energetic 
off-resonance of the excitation energy of one of the dimers 
are shown in the lower panel of Fig. [6] We flnd a general 
trend of the coupling-matrix element to decrease with 
increasing excitation energy within the data range in- 
vestigated here. A clear resonance peak of the coupling 
strength appears at the resonant coupling situation. It 
is by about a factor of 1.4 larger than the coupling at 
off-resonant situations. 



Appendix C: Off-resonant oscillation in the 
time-dependent dipole moment 



1. Electronic structure of the subsystem 

First, we discuss the electronic structure of one dimer 
and the influence of the bond length on the latter. We 
start with Na2 at the experimental bond length. For 
excitations oriented along the bond axis, Na2 is almost 
a single level atom as there is one prominent excitation 
at 2.1 eV and a second excitation at 4.1 eV with notably 
smaller oscillator strength. The changes of the excitation 



In this appendix we briefly explain how the coupling 
strength and the energetic off-resonance in a system of 
one donor molecule D and one acceptor molecule A man- 
ifest in the TD dipole moment. Note, that the resonant 
coupling case was already discussed in Ref. |HI]. Our 
considerations employ a two-state model |HZ| based on 
the assumption that the wave function of the total sys- 
tem can be separated into D and A parts due to negligible 
electronic coupling between D and A. Initially, the accep- 
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tor is in its ground state denoted by |A) and the donor is 
in an excited state |D*). This amounts to an initial prod- 
uct state |D*A) = |D*)|A) = |1). The final wave function 
corresponds to the inverse situation, where |DA*) = |2). 
Both states are characterized by their eigenvalue Ei and 
i?2- We measure possible energetic off-resonances by the 
parameter 



(CI) 



The coupling between |1) and |2) is mediated by the 
Coulomb interaction Vc- It leads to the coupling-matrix 
element 



V = (DA*|Vc|D*A). 



(C2) 



The time evolution of the two-state system with initial 
state l^'(O)) = |1) is given by 



Mt)) = aim) + ci2m) 

with the coefficients ai{t) and 02 (t) |88| . 



laaWI' = Asin^ (^^V^ + AEH^ , 



(C3) 



(C4) 



where A ~ 



and B = 



This time evo- 



lution of the coefficients corresponds to an incomplete 
oscillation with beat frequency Wbcat = yy2TAi?2 that 
depends on the coupling between the initial and the final 
state as well as the energetic off-resonance: The occupa- 
tion probability of the initial state varies around B with 
amplitude A, while the occupation probability of the final 
state oscillates with amplitude A around zero. 

The TD dipole moment di^(t) = (*(t)|f^|*(t)) of the 
acceptor, where we take the dipole operator in the 
space of the acceptor only, can be calculated as 

d^(t) = |ai(i)p(A|f^|A) + |a2WP(A*|f^|A*). (C5) 

Here, we exploited the orthogonality of \D) and \D*). 
If the static dipole moment (Ajr |j4) of A vanishes, Eq. 
(C5l simplifies to 



d^(0 = |a2(t)|'(A*|r^|A* 



(C6) 



Accordingly, with the assumption that the static donor 
dipole moment vanishes, the TD donor dipole moment 
reads 



dD(t) = |ai(i)|2(D*|fD|D* 



(C7) 



Thus, the resonance oscillation of the coefficients can be 
observed in the time evolution of the di pole mom ents 
d^(t) and d°(i). In principle both Eqs. (CgI and (|C7| 



can be used to determine the coupling-matrix element 
V . Of special importance is Eq. (C7l together with Eq. 



(C4| as a fit to the absolute of the extrema of the fc-th 
0.015 



0.01 
0.005 





donor 
fit to donor 



acceptor 
fit to acceptor 




iiilliii iiiiilii^^iiiiiliiiiiiii^iilllili 



20 40 60 80 100 120 140 
time [fs] 

FIG. 7. Donor and acceptor dipole moment (z-component) in 
a setup of two Na2 , where the bond length of the acceptor Na2 
is reduced by 0.5 bohr compared to the experimental bond 
length. We performed fits to the envelope of the oscillation 
according to the model discussed in the text. 



component of the donor dipole moment time evolution 
provides Ap^, B-pu, and Wbcat, where = (D*|f^|D*). 
It can be used to determine V ^ AE, and p^. One obtains 



V 



Apk 



Apk + Bpi 



■Wbcat 



(C8) 



and 



AE 



I Bpk 
Apk + Bpk 



(C9) 



A typical time evolution of the z-component of the ac- 
ceptor and donor dipole moments of our model system is 
depicted in Fig. [7] together with our fits to the envelopes. 
Note, that the two-level model qualitatively fits to the 
dipole oscillation of the off-resonant, coupled system of 
two molecules. However, the dipole moment envelopes do 
not perfectly follow the sin^- and cos^-shape. We under- 
stand these deviations as a consequence of the coupled 
system not being perfectly separable into D and A parts 
as we assumed in the model. Furthermore, although the 
second excitation of Na2 with polarization in z-direction 
is energetically far off and carries notably smaller oscil- 
lator strength, Na2 is not a perfect single level system, 
and therefore the system of two sodium dimers does not 
perfectly fit into the two-level model. Nevertheless, our 
approach provides a tool to determine the coupling V and 
the energetic off-resonance AE. The validity of the model 
can be checked by comparison of the thus obtained AE to 
TDDFT excitation energies of single sodium dimers with 
shifted bond length. We find good agreement and, there- 
fore, assume reasonable quality of the coupling-matrix 
element results of Fig. [6] 
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